arXiv:1507.01664vl [astro-ph.SR] 7 Jul 2015 


Draft version July 8, 2015 

Preprint typeset using D-T^;jX style emulateapj v. 08/13/06 


THE MOST IRON-DEFICIENT STARS AS THE POLLUTED POPULATION HI STARS 

Yutaka Komiya\ Takuma Suda^ and Masayuki Y. Fujimoto^^ 

Draft version July 8, 2015 

ABSTRACT 

We investigate the origin of the most iron-poor stars including SMSS J031300.36-670839.3 with 
[Fe/H] < —7.52. We compute the change of surface metallicity of stars with the accretion of inter¬ 
stellar matter (ISM) after their birth using the chemical evolution model within the framework of the 
hierarchical galaxy formation. The predicted metallicity distribution function agrees very well with 
that observed from extremely metal-poor stars. In particular, the lowest metallicity tail is well re¬ 
produced by the Population HI stars whose surfaces are polluted with metals through ISM accretion. 
This suggests that the origin of iron group elements is explained by ISM accretion for the stars with 
[Fe/H] < — 5. The present results give new insights into the nature of the most metal-poor stars and 
the search for Population HI stars with pristine abundances. 

Subject headings: stars: abundances - stars: Population HI - early universe 


1. INTRODUCTION 

The most metal-poor stars are relics of stars formed 
in the early universe. If low-mass stars with mass m < 
0.8 Mq were born from the primordial gas, they would 
still exist as nuclear burning stars in our nearby field. 
It is still controversial whether there are such low-mass, 
metal-free stars, which we refer to as Population (Pop) 
HI survivors in this paper. Recent simulations, however, 
raise the possibility that low-mass star s were formed in 
the pristine environment (see review bv iBromm 11201^ . 

Observationally, Pop HI survivors have not yet bee n 
identified despite longstanding efforts fsee lBondlll981h . 
Thanks to large-scaled surveys in the Milky Wa.y (MW) 
halo from the 1990s ([Beers &: Christlieb I l200a for re¬ 
views), several hundreds of stars with the metallicity 
of [Fe/H] < —3 are known to date and studied with 
high resolution spectroscopy (see e.g. ISuda et HI 120081 
[2ml) . Among them, six stars were found below [Fe/H] < 
—4.5 but with finite metallicity. For the most iron- 
poor star, SMSS J03I300.36-670839.3 (iKeller et al.ll20l4 
hereafter SMSS 0313-6708), iron lines are not detected 
([Fe/H] < —7.52, Bessell et al. 2015) while calcium abun¬ 
dance is determined to be [Ca/H] = —7.26 . Two other 
stars have [Fe/H] = -5.7 (HE 1327-2326. iFrebel et al.l 
l200^ and -5.4 (HE 0107-5240, iChristlieb et al.l 120021) . 
respectively. T he remaining three lie between —5 < 
[Fe/H] < —4.5 (iCaffau et al.l 120111 : iNorris et al.l l2007t 
iHansen et al.ll2014D . We refer to stars with [Fe/H] < 
—5 and —5 < [Fe/H] < —4.5 as hyper and ultra 
metal-poor (HMP and UMP) stars, respectively (c.f. 
iBeers fc Christliebl 12(3051) . 

In order to explain the absence of Pop HI sur¬ 
vivors, surface pollution with accreted metals fr om the 
interstellar matt er (ISM) has been proposed (lYoshiil 
ITmt l!i^[T^ . It is also ar gued for the extremely 
low metallicity of HMP stars ([Shigevama et al.l 120031 : 
ISnda et ^l2004 iKomiva et al.ll2009b[) . We investigated 
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the change of surface iron abundance of Pop HI survivors 
through ISM accretion using chemical evolution mod¬ 
elling wi thin the framework of hiera rchical structure for¬ 
mation ([Komiva et al.llM09bl l201fl| ). and demonstrated 
that Pop HI survivors can be polluted up to [Fe/H] ^ —5. 
The discovery of the star with metallicity smaller by two 
orders of magnitude gives us a good reason to discuss the 
range of surface metallicity of polluted Pop HI survivors. 

In this paper, we revisit the problem of surface 
pollution using an u pdated chemical evolution model 
([Komiva et al.l l20T3) . Our purpose is to explore the 
metallicity distribution function (MDF) at extremely low 
metallicity and discuss the origin of HMP/UMP stars. 
We focus on the iron abundance since lighter elements 
such as carbon and magnesium can be influenced by bi¬ 
nary mass transfer (ISuda et al.ll200l : lKomiva et al.ll2007l : 
iNishimnra et al.ll2009t) . 

2. COMPUTATION METHOD 
2.1. The Hierarchical Chemical Evolution Model 

We compute the chemical evolution of the MW within 
the framework of the hierarchical structure formation in 
the A cold dark-matter universe to evaluate the surface 
pollution of Pop HI su rvivors. The d e tailed description 
of our model is given in iKomiva et al.l ([2014] ). and we re¬ 
view the characteristics relevant to the surface pollution 
below. 

We build merger trees b ased on the ext ended 
Press-Schechter t heory ([Lacev fc Colei 119931 : 
ISomerville fc Kolat^ 119991) . The t otal mass of the 
MW is taken to be 2 x 10^^ Mq ([Li fc White! 120081 : 
iBovlan-Kolchin et al.l l2013l) . The low-mass limit of 
mini-halos is determined by the virialized temperatur e 
Tvir = 10^ K ([Tegmark et al.lfToOTt lYoshida et al.ll2003l) . 

The star formation rate is assumed to be proportional 
to the gas mass, Mgas, of mini-halos with the star forma¬ 
tion efficiency (SFE), e*, set constant. We register all the 
individual stars with metallicity Z < 0.1 Zq. We adopt 
a log-normal form of the initial mass function (IMF). 
The typical mass, mmd, of EMP stars is estimated at 
Wnid = 2.5-2OM0, much higher than the present-day 
IMF, according to the statistics of carbon-enhanced stars 
among EMP stars, which arises through the mass trans- 
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fer from AGB stars in bina ry systems (iKomiva et alJ 
[2fM[2nn^lSiida et a,l.l[2rm . We adopt rnmd = 10 Mg 
for EMP stars and 25 Mg for Pop III stars with the dis¬ 
persion of CT = 0.4, and e* = 10~^^yr~^ as fidu cial val¬ 
ues from the previous study (|Komiva et al.l[20T^ see Ta¬ 
ble [J). One of the consequences of the high-mass IMF is 
that EMP and Pop III survivors m ostly form as the sec - 
ondary members in binary systems (iKomiva et al.ll2007D . 
We take into account that the Lyman-Werner (LW) back¬ 
ground prohibit star formation in the mini-halos with 
Tvir < 10"^ K formed later than 2 = zysn = 20, al- 
though those formed earlier continue their star formation 
(iRicotti et al.ll200l: lAhn fc ShaDird[2007tl . 

We assume a uniform mixing of SN ejecta in the mini¬ 
halos. The gas outflow from the mini-halos is computed 
for each SN and HII region as a function of their ki¬ 
netic energy and the binding energy of the mini-halos. 
The conversion efficiency, rj, of the SN explosion energy 
to kinetic energy is set at 0.1. It is considered to form 
a galactic wind that enriches the intergalactic medium 
(IGM) with metals. The evolution of the metal-enriched 
wind is followed by assuming the momentum conserving 
snow-plow shell model. 

Some of the mini-halos are formed with IGM pre¬ 
enriched with metal by the galactic wind. In this study, 
we improve the treatment of the pre-enrichment by con¬ 
sidering the distance between the mini-halos, while a ran¬ 
dom spatial distribution was assumed in the previous 
work. The m-th mini-halo is regarded as pre-enriched by 
a galactic wind from the n-th halo if the distance, Dn,m-, 
between them is smaller than the radius, TZm of the metal 
enriched region around the n-th halo at formation. The 
distance Dn,m{t) is estimated following the formalism of 
the EPS theory in which the formation of a halo is de¬ 
scribed as the gravitational collapse of a spherical over- 
dense region. Two halos are considered to “merge” when 
a halo is formed by gravitational collapse with the two 
halos incorporated. In this framework, the distance be¬ 
tween the two halos that merge at time tmerge to form one 
halo of total mass Mtot may be approximated to the ra¬ 
dius, r, of the spherical region of over-density with mass 
Mtot that collapses at tmerge- 


2.2. Surface Pollution 

We trace the change in the surface iron abundance 
of each star by ISM accretion. In this work, we also 
take into account the transfer of iron from the surface 
of a primary to the secondary member in the binary 
system through the stellar wind. The accreted met¬ 
als are assumed to be mixed in the surface convection 
zones with 0.2 for giants and of 0.0035 Mg for dwarfs 
(|Fuiimoto et al.l[l995ll . 

For the ISM accretion, we adopt the Bondi-Hoyle ac¬ 
cretion formula 


rh = 


47r(Gm)^Pgas 
(K" + C.(rgas)2)3/2 


( 1 ) 


with the gas density, pgas, the sound velocity, Cs(Tgas), 
as a function of the gas temperature, and the relative 
velocity, Vr, of the star to the ambient gas, and the stellar 
mass m. 

We assume that gas in a mini-halo is concentrated in 
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Fig. 1.— The formation redshift and the mass of host mini¬ 
halos for Pop III survivors (red crosses), the second generation 
UMP stars (green squares) and the first generation UMP stars in 
the pre-enriched mini-halos (blue circles). Lines denote the mass 
pf mini-halos with the.vir.ial temperatures attached- 

its central region via isobaric contraction, i.e., 


Pgas — Pvir(^) 


-^gas 



( 2 ) 


where pvir is the averaged density of the virialized dark- 
matter halos. We set Tgas = 200 K and 50 K for 
the neutral and ion i zed p rimordial clouds, respectively 
l|Uehara fc Inutsuk^ l2000ll . When the gas metallicity 
exceeds [Z/H] = —6, we change the temperature to 
Tgas = max(10K, Tcmb(^)), where Tcmb is the temper¬ 
ature of the cosmic micr owave background, since dust 
coolin g becomes effective (lOmukai et al.l[200a ISchneideil 
1200^ . 

We set Vr = Cs before the host halo of the star un¬ 
dergoes a merger event since stars are also thought to 
be centrally concentrated. After the merger, stars from 
the smaller halo will be scattered through the merged 
halo with little dissipation of their kinetic energies, and 
we set Vr = Vcirc, where Wire is the circular velocity 
of the merged halo. The density pgas is replaced by 
Pvir(Mgas/Mh) taking into account the filling factor of 
gas clouds in the halo. 

In the case of binary systems, the accretion rate on 
each member is still poorly understood. For simplicity, 
we take an equal accretion rate for both components, i.e., 
rni/rn 2 = 1, in the fiducial model. 

The mass transfer rate in the binary systems is com¬ 
puted as a function of stellar mass and binary period with 
a given wind vel o city fro m the Bondi-Hoyle formulae as 
in IKomiva et al.l (l20C)9all . 


3. RESULTS 

Figure [T] shows the formation redshift and host halo 
mass of Pop III and UMP survivors for the model with 
the fiducial parameters given in Table [T] In this model, 
2100-2500 low-mass Pop III stars are formed in mini¬ 
halos with ~ 10® Mg around z ~ 20. Most of the second 
generation stars have metallicity of [Fe/H] > —4 since the 
iron from one SN of ^ 0.07 Mg is mixed with the gas of 
^ 10® Mg. After z = zlw = 20, newly formed mini-halos 
are prevented from forming stars by the LW radiation, 
while those formed at earlier epochs, two thirds of which 
are still metal-free, continue to form stars. The mini¬ 
halos later formed grow in mass to have Tvir > 10^ K 
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TABLE 1 

Parameters and their fiducial value 


parameter 

description 

fiducial value 


star formation efficiency 

o 

1 

1 

^md,emp 

median mass of IMF for EMP stars 

10 Mq 

^md,p3 

median mass of IMF for Pop III stars 

25 Mq 

cr 

variance of IMF 

0.4 

V 

kinetic energy fraction for SNe 

0.1 

zlw 

redshift of Lyman-Werner radiation 

20 

h 

binary fraction 

0.5 

n{q) 

mass ratio distribution of binaries 

1 

mxlrfii 

ratio of ISM accretion rate on each binary component 

1 


timefGvh 



Fig. 2.— Variations of the ISM accretion rate for the sampled 
Pop III survivors plotted against the redshift. The metallicity of 
accreting gas is color-coded. 


by z ^ 10 and form ~ 100 Pop III survivors. The sec¬ 
ond generation stars in these massive halos have lower 
metallicity of —5 < [Fe/H] < —4 because SN ejecta is 
mixed homogeneously with a larger mass of gas. The 
first generation of stars in mini-halos with pre-enriched 
IGM have metallicity around —6 < [Fe/H] < —2. They 
account for half of the UMP stars. 

Figure [2] illustrates the ISM accretion rate with the 
metallicity of accreting gas color-coded as a function of 
the redshift for the sampled Pop III survivors. The ac¬ 
cretion rate is much higher in the host mini-halos than 
in the merged halos because Vj. increases after a merger 
event. Typically Pop III survivors accrete ~ 10“® ® Mq 
of iron to have the surface metallicity of [Fe/H] ^ — 5 
for giant stars, and [Fe/H] ~ —3 for dwarfs. The sur¬ 
face metallicity of Pop HI survivors spans a wide range 
between —8 < [Fe/H] < —2 depending on the merging 
history of host halos and the mass of the primary stars of 
Pop HI survivors. The abundance pattern of the polluted 
Pop HI stars should be similar to those of EMP stars be¬ 
cause accretion in the early universe with [Fe/H] < —3 
is more efficient than metal-rich counterparts. 

Figure [3] shows the predicted MDFs for the polluted 
Pop HI and EMP stars, an d the observed samples t aken 
from the SAGA database (jSuda et al.l 1200811201111 and 
from the Hamburg/ESO (HES) survey (iSchoerck et al. I 
IM)^. The SAGA d a ta ha s more accurate metallicity 
than ISchoerck et al. I (|2009ll due to higher resolution, 
though biased toward lower metallicity of [Fe/H] < —3. 

In our previous study, we computed the MDF with¬ 
out the effects of LW feedback and SN-driven gas out- 



Fig. 3.— The predicted metallicity distribution functions 
(MDFs) compared with the observations (hisograms). Star sym¬ 
bols denote the abundances with 3D corrections for stars with 
[Fe/H] < —4.5. Left and right axes represent the observed number 
and surface number density of stars, respectively. Top panel: The 
predicted MDFs of all stars (solid line) and the polluted Pop III 
stars (dashed line) for the fiducial model (red) and models with¬ 
out SN-driven outflow (green) and without SN-driven outflow and 
Lyman-Werner feedback (blue). Bottom panel: The same as the 
top panel but for the intermediate-mass IMF model with the feed¬ 
back and the de Vaucleur profile (red) and the high-mass IMF 
model without the feedback effects (blue). Dashes and dotted lines 
denote the MDFs of Pop III giants and dwarfs, respectively. See 
text for details. 

flow (iKomiva et al.ll^IOf) . The top panel shows the de¬ 
pendences on these feedback effects. The LW feedback 
decreases the number of stars by prohibiting star for¬ 
mation in small mini-halos without metals. SN-driven 
gas outflow also decreases the number at [Fe/H] > — 2.5 
because it reduces the mass of star forming gas. For 
[Fe/H] < —2.5, however, the number of stars rather 
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Fig. 4.— The cumulative MDF of all stars (solid line) and 
polluted Pop III stars (dashed line) for the fiducial model (red) and 
the intermediate-mass IMF model (green). Histograms represent 
the cumulative nu mbers of obse r ved s tars t aken from the SAGA 
database and from IKeller et al.l 1120141') and IHansen et al.l <I2014I '1 . 
Dashed line is based on the abundances with 3D cor rections, and 
dash-d otted line is the 3D non-LTE abundance by IBessell et al.l 
II2015I 1. The model results are scaled to match the observations at 
[Fe/H] = -3. 


swells as a result of pre-enrichment. In addition, in the 
7] = 0 model, a majority of dwarf EMP stars becomes 
[Fe/H] > —3 due to efficient surface pollution. 

We compare the absolute number of stars predicted 
with the HES data considering its survey area and se¬ 
lection efficiency. The fiducial model predicts a smaller 
number of EMP survivors, while the model without the 
feedbacks gives the number in agreement with observa¬ 
tions when we assume a homogeneous surface density 
of EMP halo stars. Additionally, if we adopt the de 
Vaucouleur density profile for the stellar halo, the scal¬ 
ing factor can be smaller by a factor of five. It is to 
be noted, however, that the number of EMP survivors 
is determined by the production ratio of low- to high- 
mass stars since the latter dominates the chemical evo¬ 
lution. In our models, it depends on the nimd and a of 
the IMF, the binary fraction ft, and the MRD. For ex¬ 
ample, we present the result with the intermediate-mass 
IMF, TOmd,emp = 3 Mq and TOmd,p 3 = SMq, ill the bot¬ 
tom panel. This model results in the number of survivors 
larger by 0.7 dex, compatible with the observation un¬ 
der the de Vaucouleur density profile. The shape of the 
MDF has little dependence on the adopted IMF. In all 
cases, stars with [Fe/H] < —5 are polluted Pop HI stars 
except for a small contribution from the first generation 
stars in pre-enriched halos. 

Figure [4] shows the cumulative MDFs. The ob¬ 
served dat a are taken from the SAGA sample, and 
added from iKeller et al.l (|2014D : IHansen et aH (|2014D and 
IBessell et al.l ( 20I5I) to complete the HMP/UMP sample. 
We also show the abundances with 3D corrections when 
available. The broad distribution of polluted Pop HI 
survivors agrees with the distribution of the observed 
stars at the lowest metallicity range. The metallicity 
range of the HMP stars including SMSS 0313-6708 is 
covered by the polluted Pop HI survivors. We estimate 
the iron abundance of SMSS 0313-6708 to be around 
[Fe/H] ~ —7.58 from the observed calcium abundance 
([Ca/H] = —7.26) and the typical enhancement of cal¬ 
cium for EMP stars ([Ca/Fe] ~ -1-0.32). The iron abun¬ 


dance is unlikely to be [Fe/H] < —8.0 in our scenario 
because 98% of EMP stars show [Ca/Fe] < -1-0.74. Note 
that the expected number of observed stars is ^ 1 at 
[Fe/H] < -7.0, - 0.6 at [Fe/H] < -7.6, and - 0.4 at 
[Fe/H] < —8.0. We see a slight overabundance for the 
model results around [Fe/H] ~ —5 to —4. This is an arti¬ 
fact resulting from the assumption of homogeneous mix¬ 
ing of SN ejecta in such massive halos as Mh > 10^ Mq, 
and can be alleviated by considering possible inhomoge¬ 
neous mixing. 

We also check the dependencies on other model param¬ 
eters. There are minor corrections to the fiducial model 
by changing rhxjrhi. Binary parameters such as ft and 
n(g) affect the number of Pop HI and EMP survivors but 
have little effect on the metallicity distribution. On the 
other hand, the MDF at the lowest metallicity is affected 
by the SFE; we find a rapid decline at [Fe/H] < —6 for 
e* = 10“^°yr“^, because of the higher metallicity that 
the ISM reached before the merger events. 

4. CONCLUSIONS AND DISCUSSION 

We have explored the change in the surface iron abun¬ 
dance of low-mass survivors of Pop HI stars by the accre¬ 
tion of ISM using the chemical evolution model within 
the framework of hierarchical structure formation. 

Our model can reproduce the MDF observed for Galac¬ 
tic halo stars with [Fe/H] < —2.5. The MDF of polluted 
Pop HI giants peaks around [Fe/H] ~ — 5 and extends 
to cover the whole metallicity range of HMP stars in¬ 
cluding SMSS 0313-6708. A wide spread of metallicity 
of polluted Pop HI stars arises from the intervals from 
the birth to the merger of host halos and the mass of 
primary stars. The results little depend on the adopted 
IMF, and the binary parameters, but are dependent on 
the SFE. 

Surface pollution is most effective in the mini-halos 
of small mass because the Bondi-Hoyle accretion rate 
is strongly dependent o n the relative velocity of stars 
to the ambient gas fe.g.. lSuda et al.l[2004 iKomiva et al. 
l2009al) . This explains the difference from iFrebel et ah 
(|2009fl . who considered only the s urface pollution dur¬ 
ing pa ssage through the MW disk. I.Iohnson fc Khochfarl 
(l2ml) argued the suppression of ISM accretion by al¬ 
leged solar-like winds from Pop HI survivors. However, 
little is known about the mass loss from Pop HI stars 
with a complete lack of metals and also in the possible ab¬ 
sence of magnetic fields. In addition, their argument on 
the dynamics of Pop HI survivors in the mini-halos may 
not be relevant for those bel onging to binaries wit h more 
massive stars. We note that iHattori et al.l (|2014D report 
possible observational evidence for the surface metal pol¬ 
lution, finding that the median metallicity of halo G-type 
dwarfs may be systematically higher than that of K-type 
dwarfs with the same kinematics. 

The present results indicate that such small abun¬ 
dances of iron group elements as observed from HMP 
stars are readily explained in terms of the ISM accre¬ 
tion. It is true that HMP/UMP stars show various pe¬ 
culiar abundance patterns other than the large enhance¬ 
ment of carbon, i.e., the enhancement of N, O, and light 
elements (Na, Mg, and Al), and also of Sr. These el¬ 
ements can be synthesized and brought to the surface 
in low- and intermediate-mass stars during the AGB 
phase ([Nishimura et al.ll2009t lYamada et al.l 120151) . and 
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transferred onto t heir low-mass secon dary companions 
(jSuda et al.ll2004 iKomiva et alJl2007ti . We will discuss 
the relevance of the nucleosynthesis and material mixing 
in AGB stars to t he observed abun dance patterns in a 
subsequent paper (iSuda et al.ll2015[l . 

As an alternative scenario for the origin of HMP 
stars, it has been argued that they are the stars formed 
from the matter polluted with the ejecta from “faint 
SNe” with large C/Fe r a tios le.g. lUmeda fc Nomotol 
Iwamoto et 'al]l2005ll . iKeller et al.l (I 20 TII : llshigakil 


( 201^ : lMarassi et al.l '( 201^ . and iBessell et alTlI^Olsil ar- 


gue that the abundance pattern of SMSS 0313-6708 is 
consistent with the ejecta of a peculiar SN with small 
(or no) iron ejection. However, this scenario demands 
manipulation of reducing the ejection of iron by many 
orders of magnitude 10^-th or less). It must be noted 
that large C/Fe ratio results solely from the reduction 
of iron ejecta in “faint SNe” since the amount of car¬ 
bon produced hardly differ from ordinary SNe. There 
should be an account for producing such large carbon 
abundances and variations as [C/H] = —1.3- — 2.4 ob¬ 


served for H MP stars with almost similar amount of car¬ 
bon ejected (jSuda et al.ll2004 iKomiva et al.ll2007|l . 

In recent simulations of first star formation, low-mass 
stars can be formed at the same sites as massive stars 


('e.e. IClark et al.l 120081 

Greif et al.l 120111: ISmith et alJ 

120111 IDoncke et all 201.8 

: Stacv & Bromml 120141). This 


opens a new path for low-mass Pop HI star formation 
in addition to ioniz ed primordial gas in the mini-halos 
with Tvir > 10^ K (|Ueliara fc Imitsui^l200flf l. Our sce¬ 
nario indicates that Pop HI survivors, found in the Milky 
Way halo, should have suffered more or less from surface 
pollution by the ISM accretion. In order to confirm the 
formation of Pop HI low-mass stars, therefore, we have 
to seek for un-polluted Pop HI survivors, which are ex¬ 
pected to have been expelled from the host halos. We 
will discuss the search for these pristine Pop HI stars in 
a forthcoming paper. 


This work is partially supported by Grant-in-Aid for 
Scientific Research (23224004, 25400233, 25800115). 
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